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We report theoretical and experimental investigations of the flow instability responsible 
for the mushy-layer mode of convection and the formation of chimneys, drainage channels 
devoid of solid, during steady-state solidification of aqueous ammonium chloride. Under 
certain growth conditions a state of steady mushy-layer growth with no flow is unstable 
to the onset of convection, resulting in the formation of chimneys. We present regime di- 
agrams to quantify the state of the flow as a function of the initial liquid concentration, 
the porous-medium Rayleigh number, and the sample width. For a given liquid con- 
centration, increasing both the porous-medium Rayleigh number and the sample width 
caused the system to change from a stable state of no flow to a different state with the 
formation of chimneys. Decreasing the concentration ratio destabilized the system and 
promoted the formation of chimneys. As the initial liquid concentration increased, onset 
of convection and formation of chimneys occurred at larger values of the porous-medium 
Rayleigh number, but the critical cell widths for chimney formation are far less sensitive 
to the liquid concentration. At the highest liquid concentration, the mushy-layer mode of 
convection did not occur in the experiment. The formation of multiple chimneys and the 
morphological transitions between these states are discussed. The experimental results 
are interpreted in terms of a previous theoretical analysis of finite amplitude convection 
with chimneys, with a single value of the mushy-layer permeability consistent with the 
liquid concentrations considered in this study. 



1. Introduction 

The solidification of binary alloys and the formation of so-called mushy layers are 
importa nt in applications ac ross geophysics, geology and industry, ra nging from sea ice 
grow th ( Feltham et all [20061 ). to solidific ation in the Earth's interior ( Bergman fc Fearnl 



1994 ). to the formation of metal castings ( Copley et aLlll970h . Mushy layers are partially 



solidified reactive porous materials comprised of solid crystals and interstitial melt, and 
formed as a r esult of morphological in stability during the solidification of a two com- 
ponent melt ( Mullins fc Sekerka 1964) . An issue of significance in many applications 



is understanding the conditions for both the occurrence a nd the resultin g dynamics of 



buoyancy-driven convection within the mushy region (see IWorsterl 119971 for a review.) 
Convective flow in a mushy layer transports fluid into regions of differing concentra- 
tion, which can lead to local solidification or dissolution of the solid matrix and changes 
in the effective material properties. The permeability of the porous matrix depends on 
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the local solid fraction, and hence flow focussing can accelerate the growth of the un- 
derlying convective-flow instability, and eventually lead to the formation of chimneys - 
channels of zero solid fraction that provide efficient pathways for buoyancy-driven flow 
out of the mushy layer. Understanding the onset and dynamics of convection in mushy 
laye rs with chimneys is of importance in geop hysics for predicting brine fluxes from sea 
ice ( Wettlaufer et al. 19976 ; Wells et al. 201 ll ). and in metallurgy because material het 



eroge neities associated with chimneys can lead to defects in metal castings (jCoplev et al 
197(f) . 



Experimental studies of chimney formation in mushy layer growth typically take one of 
two forms. Numerous so-called fixed-chill experiments have studied the transient growth 
of a mushy layer from a cold isothermal boundary and show that ch imneys typically 



form after the mushy layer exceeds a critical thickness (see, for example, lAussillous et al 



20061 and references therein). As the thickness of the mushy layer continues to increase, 
it is observed that the mean sp acing between chimneys also increases, with extin ction 
of the flow in certain chimneys ( Wettlaufer et aZ.1 Il997al ; ISolomon fc Hartley 1998 ). Al- 
ternatively one also considers steady-state mushy layer growth. Directional solidifica- 
tion experiments, where the sample material is translated between two heat exchangers, 
indicate that multiple chimneys form for growth r ates slower than a critical thresh- 
old that depends on the initial liquid concentration ( Peppin et all [2008 ; Whiteoak et al 



20081) . Neufeld fc Wettlauferl (|2008af ) used growth from a boundary with controlled time 
dependent temperature in order to achieve steady-state growth and study the influence 
of an external shear flow on the onset of chimney formation. For both transient and 
steady state growth the critical threshold for chimney convection can be interpreted in 
terms of an appropriate mushy-layer Rayleigh number exceeding a critical value, so that 
buoyancy effects overcome dissipation and lead to the onset of convective fluid flow. 

Theoretical modelling of convective flow in mushy layers has been approached from a 
variety of angles. Linear and weakly nonlinear stability analyses have considered the on- 
set of convect ion during directional solidification under a wide variety of conditions (see, 
for example, IWorsterl 1 19971 iGuba fc Worsterl l2010l : Roper et al. 2011 . and references 
therein for a review and more recent summaries.) Several studies have revealed that 
the bifurcation to both steady and oscillatory convective flow states can be subcriti- 
cal ([Amberg fe Homsvlll993t lAnderson fe Worsterlll995t IGuba fe Worsterll2006l) depend- 
ing on the growth conditions. Hence, this raises the possibility that finite amplitude 
perturbations may allow states with chimney convection to be accessed even when the 
base state of no flow is linearly stable. Convective states of finite ampl itude have been con- 
sidered either in direct simulations of particular growth conditions (IFelicelli et aZ.M l998 



Bcck ermann et al. IbOOOllGuo fe Beckermannl200alHeinrich fe Poirierl2004HOertling fe Wattsl 
2004t Ijain et a/.ll2007l : iKatz fe Wor stcr 2008 1) or via studies of the system dynamics with 
an assumed perio dic array of chimneys (jSchulze fe Worster | ll998HChung fc Worsterll2002t 



Wells et allmmh . In a study of particular relevance here IWells et al\ (|2010l ) used a nu- 
merical approach to consider the influence of changes in the spacing between chimneys on 
the evolution of system dynamics and the stability of solutions with steady-state chim- 
ney convection. This approach demonstrates that the observed convective flow states are 
suppressed via a saddle-node bifurcation as the chimney spacing is reduced. 

In this study we focus on the influence of a finite sample width on the formation of 
chimneys during mushy layer convection. I n w e describe how the theoretical model 
of directional solidification of IWells et al\ (|2010h can be used to study the dynamics 
and stability of nonlinear convective flow states of finite amplitude. An experimental 
procedure is described in Sj3]that pr oduces steady-state mushy layer growth for compari- 
son to the theoretical predictions of Wells et all ( 2010h . An extensive set of experiments 
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Figure 1. Schematic of the theoretical framework. Liquid of initial concentration C and tem- 
perature Too is translated downward at speed V between hot and cold heat exchangers, and 
partially solidifies to form a mushy layer of thickness h(x,t). We adopt a local co-ordinate sys- 
tem (x,z) so that a eutectic solid of temperature Te and concentration Ce forms at z = 0. We 
solve in a domain < x < I and < z < h(x,t) as indicated by the dashed outline, with a 
chimney of width a(z,t) located at x = located in this domain. With suitable boundary con- 
ditions, described in the main text, this domain describes either a periodic array of convection 
cells separated by chimneys of width 2a(z, t) centred at x — and with inter-chimney spacing 
21, or else a single convection cell with chimney of width a(z, t) against an impermeable free-slip 
vertical boundary. 



described in 21 allows us to identify conditions for chimney formation for a variety of 
sample widths, growth rates and liquid concentrations, with comparison to the theoreti- 
cal results providing an indirect estimate of mushy layer permeability for several different 
liquid concentrations. We conclude with a discussion in SJSJ 



2. A model of mushy layer convection with chimneys 

A t heoretical model of mushy layer convection with chimneys was used bv I Wells et 



( 2010f ) to investigate the dynamics of nonlinear convective flow, and is applied here 



to determine conditions under which flow states with a chimney can exist. We briefly 
summarise the essential features of the model below, noting that further deta i ls of t he 
governing equations and boundary conditions are elaborated on in IWells et ali (l2010lh 



2.1. Physical and theoretical framework 

The model considers the directional solidification in two dimensions of a binary solution 
of local temperature T and local concentration C, translated at a speed V between hot 
and cold heat baths, as illustrated in figure [TJ This results in a porous mushy region of 
thickness h(x, t) bounded by an impermeable eutectic solid with temperature Te and 
and concentration Ce at z = 0, and an overlying fluid for z > h with assumed uni- 
form composition C and far- field temperature T^. We simulate flow with a chimney 
and single convection cell of width I, utilizing boundary conditions that simultaneously 
allow us to consider both convection cells that are either isolated, or in periodic arrays. 
The temperature T, local solid fraction <j) and Darcy velocity u wi thin the mushy re- 



gion are calculated using time-dependent ideal mushy layer theory (jSchulze fc Worster 



20051) which conserves heat and salt, accompanied by a condition of local thermodynamic 



equilibrium which couples the temperature and liquid concentration via a freezing point 
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relationship T = Tl(C). We use a linear approximation Tl(C) = Te + T(C — Ce) for the 
freezing point relationship, where V is constant. Fluid velocities are determined from the 
incompressible Darcy flow law, where the permeability is assumed to be II = IIo(l — 4>) 3 
for some reference permeability IIo. We apply symmetry conditions on T and no normal 
flow at x = I, equivalent to either an insulating, impermeable boundary or a periodic 
array of convection cells. 

The pure liquid flow in the narrow chimney is driven by b uoyancy forces or pressure 
gradients and is parameterised using lubrication theory (see IWells et all 2010l . and ref- 



erences therein) where we apply symmetry conditions on the velocity u at x = 0. This 
velocity boundary condition can describe either flow within two convection cells sepa- 
rated by a chimney of width 2a(z, t) centred at x = 0, or else a single convection cell with 
chimney of width a(z,t) against an impermeable free-slip vertical boundary (a discussion 
of the likely influence of using a free-slip condition rather than a no-slip condition follows 
at the end of i )2.1[) . Combining the lubrication approximation with similar asymptotic 
approximations for heat transport yields nonlinear boundary conditions on the left hand 
side of the mushy region, which is treated as a singular interface at x — 0. These con- 
ditions depend on the width of the chimney, with a(z,t) updated to satisfy a condition 
of marginal equilibrium to just remove any constitutional supercooling. The boundary 
conditions at the mush liquid i nterfac e and procedure for updating h(x,t) are described 
in more detail in Wells et al. ( 2010( ). A pertinent point is that this approach uses a 



boundary layer approximation that describes the circulation and heat transport in the 
overlying fluid that is induced by the mushy-layer mode of convection, but neglects any 
component of buoyancy- driven convection that arises due to density differences within 
the overlying fluid region. 

This model is characterized by six dimensionlcss parameters, 

.VI PogPjCo - C E )Tlo _ Cs-Cp 

k /j,v Go — Ce 

9 _ T oo -T L (c ) l T,,-IbYl ran 

where the fluid density is approximated by p$ [1 + j3(C — Ce)] for constant haline coeffi- 
cient P and reference density po , /i is the dynamic viscosity, g is the gravitational acceler- 
ation, k is the thermal diffusivity, c p is the specific heat capacity, L is the latent heat of 
fusion, and we have made the simplifying assumption that material pr operties are uniform 



in bo th solid and liquid phases in line with ideal mushy layer theory (jSchulze &: Worster 



20051) . The numerical results presented herein investigate how the dynamics change under 
variations of the dimensionless sample width A and mush Raylcigh number R m which 
characterizes the strength of buoyancy effects compared to dissipation. The experimental 
concentration and temperature conditions are reflected in the concentration ratio 9oo, 
and 5, with the Darcy number Da representing the relative permeability of the porous 
matrix. It is difficult to directly constrain and measure the permeability in experimental 
settings, and so in ij4.2l we determine a reference value of IIo = 1-6 x 10~ 5 cm 2 for a 
range of liquid concentrations, by comparison between experimental data and theoretical 
prediction. This provides an indirect estimate of the permeability for solidification in 
a Hclc-Shaw cell. Comparison to other estimates in different settings will be discussed 
in Jj5] Importantly, we note that an independent set of simulations that vary Da with 
all other parameters held fixed show that whilst th e chimney width va r ies wi th scaling 
aV/n ~ Da 1 / 3 consistent with the scaling analysis of Schulze fc Worster ( 1998|) . the fluid 



flux remains constant at leading order over the range 10 4 < Da < 5 x 10 3 . The same 
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features have also been observed in a recently-developed simplified model of convection 
in mushy layers (Rees- Jones & Worster, personal communication). This behaviour is 
consistent with the flow dynamics being controlled by buoyancy-driven flow within the 
mushy region, whilst the chimney width adjusts passively to accommodate the necessary 
liquid flux. Hence, the mushy-region flow dynamics do not depend strongly on the ratio 
of flow resistance in the porous mushy region to viscous resistance in the liquid region, 
consistent with only a weak dependence on Da. We maintain the value Da = 5 x 10~ 3 
throughout all numerical results reported here. 

The fact that flow in the mushy region is relatively insensitive to viscous resistance 
in the liquid region also suggests that qualitatively similar results would be obtained by 
application of a free-slip and a no-slip boundary condition in the chimney. The boundary 
condition describing flow at the chimney edge takes the form 



Da 



aV 



dip 



) 3 dx 
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013a 



chp_ 
dx ' 



(2.2) 



3R m (l 

dWells et ffll.ll2010h . where the streamfunction ip satisfies u = ipx) and the constant 

pre-factors are ct\ = 011 = 013 = 1 for a free-slip boundary, or a\ = 1/4, CX2 = 2/3 
and a-i = 1/2 for a no-slip boundary. For fully developed chimneys, with aV/n ~ Da 1 ^ 3 
the fin al term on the right-han d side of (|2.2p is 0{Da}^) smaller than the remaining 
terms (jSchulze fc Worsted If 9981 ). and so an order-one change in 03 does not impact the 
leading order behaviour. Next, we note that the numerical coefficient of the buoyancy 
term (featuri ng q 2 ) has been d etermined using a Polhausen approximation in the current 
calculations ( Wells et gl. 201 Oh . Previous calculations using a different form of approxi- 
mation yie ld a different numerical c oefficient, but qualitatively similar mushy-region-flow 
dynamics ([Chung fc Worsterll2002l ). Hence, we expect an order-one change in a 2 to also 
produce qualitatively similar dynamics. Finally, the difference in ct\ between free-slip and 
no-slip boundaries can be compensated for by a rescaling of the Darcy number Da — > 4Da. 
Hence, because the flow in the mushy layer is insensitive to Da at leading order, we expect 
no qualitative difference in the leading-order dynamics within the mushy region between 
cases with a free-slip or a no-slip boundary condition in the chimney. 



2.2. Numerical method 



The solution procedure combines two different numerical approaches involving both time 
dependent calculations, and an arc-length-continuation scheme. The time dependent 
simulations combine second-order finite differences in space, with semi-imp licit Crank 



Nicho lson timcstepping, and free boundaries updated via relaxation (see IWells et 



2010L for further details). The result ing system s for temperature and streamfu nction 
were solved using a multigrid iterative method ( Adams 19891 : Briggs et al. 2000l). This 
time dependent code was complemented by an arc-length-continuation method ((Keller 
19771) that traces both stable and unstable solution branches, allowing us to verify that 
the relaxation procedures do not have a strong influence on the stability of the system. 
In order to reduce the computational costs in solving the systems involving Jacobian 
terms that are required for arc-length continuation, the predictor-corrector calculations 
are computed from projections of the finite difference solutions onto a basis of products 
of Chebyshev polynomials in x and y. This allows the main characteristics of the solution 
to be characterised by a smaller number of degrees of freedom. The resulting solutions 
predicted by arc-length continuation in Chebyshev space are then refined using a steady 
version of the finite-difference code. 
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Figure 2. Variation of the calculated dimensionless brine flux F with convection cell width A 
for steady state flow with R m = 32.5 "W = 11.7, S = 4.8, 6^ = 0.185, and Da = 0.005, consis- 
tent with solidification of a 27wt% aqueous ammonium chloride solution described in ij3] The 
system shows hysteresis, with multiple stable steady states. The stable upper branch of chimney 
convection (red circles) and stable lower branch of no flow (blue squares), are accompanied by 
an intermediate unstable branch of chimney convection (black crosses). The system undergoes a 
saddle-node bifurcation at (A c , F c ), with the stable state of chimney convection ceasing to exist 
for A < A c . Failure of convergence of the iterative corrector scheme prevents further continuation 
along the unstable branch. 



2.3. Stability boundaries 

In this paper we are interested in the influence of sample width on the formation of 
chimneys in a mushy layer. To determine conditions for the existence of states with 
chimneys, the arc-length-continuation method is used to trace solution br anches as A 
varies with R m , ^ , S, 9qo and Da all held fixed, applying the method of IWells et all 
to a region of parameter space relevant to our experiments. Figure [2] shows a 
typical example of the system evolution, where F is the dimensionless solute flux from 
the chimney. Two steady states are observed, an upper branch of chimney convection (red 
circles) and a lower branch of no flow (blue squares) , along with an intermediate unstable 
branch of chimney convection (black crosses). Starting on the upper branch of stable 
chimney convection, as A is decreased the system undergoes a saddle-node bifurcation at 
A = A c so that the chimney convection state ceases to exist for A < A c . Having determined 
the critical wavelength A c , the above procedure is repeated with varying R m but 5, 
0oo and Da held fixed to obtain a stability boundary of the form A = A c (R m ). Examples 
of such stability boundaries and their comparison to experimental data are shown later 
in figure [7] and discussed in f]4] Note that, by construction, the theory predicts that flow 
with a single chimney at the centre of two symmetric convection cells requires twice the 
domain width as the corresponding state with a single convection cell and chimney at 
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Figure 3. (a) Schematic of the experiment apparatus. See text in Sec. 3.1 for detailed identifi- 
cation of the various components, (b) Phase diagram of the NH4CI-H2O solution. The vertical 
dashed lines indicate the initial liquid concentrations chosen in the experiment. Solution den- 
sity isopleths are indicated b y the faint lin es with positive slope above the liquidus and were 
constructed using data from l|Councilll2003f ) 



the sidcwall. Hence one intuitively expects that a state of convection with a single wall 
chimney can be achieved at smaller Raylcigh numbers and narrower cell widths than 
convection with chimneys away from the sample walls. 



3. Experimental apparatus and methods 

3.1. Experimental apparatus 

A schematic diagram of the apparatus is shown in Fig. [3^. A Hclc-Shaw sample container 
(A) of inner dimensions 25.40 x 15.24 x 0.32 cm is surrounded by an air-flow vessel (B). 
The container walls (C) are made of 0.8-cm thick, optically clear plexiglas, and extend 0.6 
cm below the top surface of the bottom plate (D). The seal which confines the working 
fluid is provided by high-strength silicon adhesive (RTV 732) below the walls. Spacers 
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(E) made by low-thermal-conductivity plexiglas divide the container into a subset of 
cells with adjustable width I. The width of the spacers l s =1.3 cm is chosen to provide 
sufficient thermal insulation between the convection cells. The thermal diffusion time 
over the spacers, l 2 s /n s is of order hours (k s is the thermal diffusivity of plexiglas), 
and is comparable to the whole time span of the solidification process. In this manner, 
mushy-layers grow in each of the convection cells independently, and each is subjected 
to spatially uniform cooling from the bottom plate. 

The cell bottom (D) is made from a 1.3-cm thick high thermal conductivity copper 
plate. Three thermistors are embedded in (D) with a 0.3-cm distance to the top plate 
surface (shown by the solid black circles in Fig. [3^,). These monitor the bottom temper- 
ature of the sample Tb- To inhibit chemical reactions with the working fluid, the top 
surface of (D) is painted with an anti-corrosion coating (PRO 15) 0.1mm thick. When 
the highest thermal flux is applied through the sample, this coating is responsible for a 
temperature difference of 0.3K between Tb and the fluid temperature measured at the 
bottom of the sample. Corrections have thus been made for Tb to accurately reflect the 
bottom fluid temperature. 

The sample cell is enclosed by an air-flow vessel (B). Ambient air enters and leaves 
through outlets in this vessel, as indicated by the arrow heads in Fig. [3^,, at a rate of 
approximately 50 cm 3 /s. The air flow has been thermally stabilized at = 20 ± 0.2°C 
in a large reservoir before entering the vessel to provide a steady thermal environment. 
Simultaneously we flow dry air flow through this region to suppress water condensation 
on the sample walls (C) during the solidification process. An aluminium top plate (F) 
maintains a constant far-field temperature Too of the sample. The working fluid confined 
in the cell is an aqueous NH4CI solution with initial concentration Co. The density and 
temperature of the solution are measured before the experiments using a oscillating-type 
densitometer (Antor Paar, DMA 35N). The concentration Co is then chosen between 
23wt% and 27wt% with an accuracy better than ± 0.1wt% using the isopycnal lines as 
shown in the phase diagram in Fig. \3jp. 

In thin cells the thermal boundary conditions of the sample may be affected by heat 
flow through the walls to t he working fluid, w hich has been considered explicitly in a 
recent mathematical model (jPeppin et aZ.ll2007l ). Because of the finite thermal resistance 
of the wall material, a sample cell with wall temperature profiles that match entirely 
the fluid temperature is necessary to eliminate the lateral heat flow effects. The double- 
sided wall construction in the present experimental system effectively minimizes this heat 
loss effect. For a quantitative evaluation of the lateral heat lost through the sidewalls, 
we measure the fluid temperatures by moving a temperature probe along the shortest 
dimension of the cell, at a horizontal level 3mm above the mush-liquid interface before 
the formation of the chimneys. The temperature variation is within 0.03K over the full 
thickness (0.32mm) of the cell, with the centre fluid temperature slightly higher than 
that near the sidewalls. Thus, wc conclude that the vertical heat flux through the sample 
dominates the lateral heat lost along the sidewalls, and that the temperature field in the 
liquid region is essentially two-dimensional. 

Underneath the bottom of the sample cell there is an array of Peltier devices (G) (ther- 
moelectric coolers) sandwiched between the copper plate (D) and a 2.5-cm thick copper 
manifold (H). The temperature of the manifold is regulated through the circulation of 
anti-freezing coolant. High-conductivity thermal compound is pasted between the Peltier 
devices and the copper plates to enhance the thermal contacts and improve the thermal 
homogeneity on the bottom plate of the sample. Under typical experimental conditions, 
the horizontal temperature heterogeneity along the sample bottom is within 0.05K, as 
measured by the three thermistors embedded in part (D). 
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Supporting the copper manifold (H), a thermal insulating based shield (I) rests on a 
levelling plate (J) through which the entire apparatus can be levelled within 0.005 radian. 

3.2. Thermometry and thermal control protocol 

The bottom-plate temperature Tb and top-plate temperature T x are measured with ther- 
mistors imbedded in the bottom plate (D) and the top plate (F). The fluid temperature 
in the liquid region is measured by thermistors inserted inside the sample through 0.6- 
mm diameter stainless steel tubing. All thermistors are calibrated to a precision of lmK 
against a Hart Scientific model 5626 platinum resistance thermometer and are traceable 
to the ITS-90 temperat ure scale. Consistency be tween the thermistors is within 0.01K. 
Recent experiments (jPeppin et all 12007 , 120081) investigated steady state mushy-layer 



growth using an experimental apparatus in which the solution is translated at prescribed 
rates between two heat exchangers. In the present experiment, steady growth is accom- 
plished by time dependent temperature control as Tg = T B — jt where 7 is the constant 
cooling rate. 

Following the protocol developed in a recent experimental study, c ontrol of the sample 
bottom -plate temperature Tb is achieved by two cooling mechanisms (jNeufeld fc Wettlaufer 



2008q( ). First, we use an array of Peltier devices (part (G) in Fig. [3^) to produce a ther- 
mal flux from the bottom plate (D). The downward heat flux, which is proportional 
to the applied electrical power to the Peltier devices, is controlled in a feedback loop 
that incorporates the instantaneous temperature readings for Tb , and is delivered to the 
Peltier devices by a programmable power supply (Sorenson, Model DCS 40-25). A sec- 
ond ingredient of the cooling-control settings is the coolant-circulating copper manifold 
(H). Anti- freeze coolant is circulated in a double-spiral channel built inside manifold (H), 
driven by two Neslab (RTE-7) baths at a flow rate of 500 cm 3 /s and a cooling capacity 
of 1000 Watt at 20° C. The bath-temperature stability achieved over long time scales is 
within 0.01°C, as measured by the thermistor embedded in the manifold (H). The Peltier 
devices provide fast control of the temperature variation on the bottom plate (D). Dur- 
ing the experiment a prescribed linearly decreasing bottom-plate temperature Tb of the 
sample is maintained down to approximately -40°C (Fig.|4j:), which forces steady growth 
of the overlying mushy layer as depicted in Fig. HJd. 



4. Results 

4.1. Steady state mushy-layer growth 

One representative set of data for the mushy-layer growth as a function of time is shown 
in Fig. [4)3. There are three distinct periods, with different characteristic variations of h(t) 
with time. To the left of the first dashed line at Tb =T\ = — 17.6°C (shown in Fig. Hb), 
the thickness of the mushy layer increases linearly with a growth speed of 6.60/im/s. No 
eutcctic layer is formed in this period. Formation of a eutectic solid layer underlying the 
mushy-layer starts at Tb = Ti. As the bottom-plate temperature Tb further decreases 
an initial transient-growth period appears between (— 17.6°C, — 24.4°C), as indicated by 
the two vertical dashed lines. Growth of the mushy-layer in this period is dependent on 
initial conditions, particularly the cooling history before Tb reaches T\. 

To the right of the second dashed line at Tb = T\ = — 24.4°C, both the advancing 
speed of the mush-liquid interface, Hi, and the eutectic- mush interface, He, are constant 
and equal V = 4.86/im/s. The thickness of the mushy-layer, h = Hi — He, thus remains 
a constant independent of the initial cooling process. The thickness of both the mushy 
layer and the eutectic layer in this steady-growth period are shown in Fig. [4^. 

From the description above it is clear that the advancing speed of the mush-liquid 
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Figure 4. (a) Image of a mushy-layer growing from a 27wt% NH4CI solution with cooling 
rate 7=0.5K/min. The bright white curve is the eutectic-mush interface. The horizontal black 
bar is 1 cm. (b) Time-recording of the mushy-layer thickness depicting three growth periods. 
Black curve: height Hi of the mushy-liquid interface; Blue curve: height He of the eutectic 
solid layer; Red curve: thickness of the mushy-layer h — Hj — He- (I) Linear-growth period. 
In this period the bottom-plate temperature Tb is above the eutectic temperature Te- The 
mushy-layer grows linearly as a function of time until the eutectic solid forms. (II) The initial 
transient-growth period (Ti > Tb > Tb). The eutectic layer starts to grow in this period. The 
mushy-layer thickness depends on time and finally converges to a constant. (Ill) Steady-growth 
period {Tb < T2). The mushy-layer thickness remains constant, (c) The linearly decreasing 
bottom- plate temperature Tb- The two dashed lines in (b) and (c) show the two temperatures 
at state transitions: Ti = -17.6°C, T 2 = -24.4°C. 



interface V is not a prescribed control parameter in this study but depends on the 
bottom-plate cooling rate 7 and the initial liquid concentration Cq. In Fig. [5^ we show 
that at high cooling rates, when 7 > l.OK/min, V increases with Co, whereas at low 
cooling rates 7 < l.OK/min, V is nearly independent of Co- 
in Fig. [5b we show data for the mushy-layer thickness h during a steady state of 
chimney convection at cooling rates covering 0.5K/min < 7 < 2.5K/min. The dependence 
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Figure 5. (a) Growth velocity V of the mush- liquid interface as a function of the cooling rate 7 in 
a steady-growth state with / = 2.54 cm. Black squares: Co = 27wt%; Red circles: Co = 26wt%; 
Green diamonds Co = 25wt%; Blue triangles: Co = 24wt%; Dark yellow triangles: Co=23wt%. 
Dashed lines: spline fitting curves for each set of data points, (b) Height of the mushy-layer ftasa 
function of the cooling rate 7 in a steady-growth state with Co = 27wt%, I — 2.54cm. The error 
bars show the variation in h at different runs under the same experimental conditions. Inset: 
The mushy-layer height h as a function of Tb at three cooling rates. Red curve: 7=0.5K/min; 
Black curve: 7=1.0K/min; Blue curve: 7=2.5K/min. 



of h as a function of the bottom-plate temperature is shown in the inset for three cooling 
rates 7 =0.5, 1.0 and 2.5K/min. When Tb is sufficiently low the system easily reaches 
a steady state. The mushy-layer thickness h is then evaluated when it converges to a 
constant and is independent of Tb- Figure [SJj shows that the mushy- layer tends to be 
thicker at lower 7, leading to a higher mu shy- layer Ray le igh number. These trends agree 
qualitatively with theoretical predictions (|Worsterlll992l) . 

Results for h at the same cooling rates 7 are repeatable from run to run and are 
independent of the pre-cooling process. The inter-experimental variation of the steady- 
state mushy-layer thickness h is within 5%, as indicated by the error bars in Fig. [SJd. 
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Experiments with very low cooling rate 7 <0.5K/min show that the advancing speed 
of the eutectic-mush interface is always larger than the mushy-layer interface (He > 
Hi), leading to a decreasing mushy-layer thickness and the system does not reach a 
steady state within the experimental timeframe. Due to the finite cooling capacity of 
the experimental apparatus, at experiments with high cooling rate 7 >2.5K/min, the 
bottom-plate temperature Tg deviates from the usual linearly decreasing trends in the 
low-temperature range. Data analysis in these experiments is limited by the shorter 
period of growth. 

4.2. Instability and onset of mushy-layer mode convection 

For mushy-layer Rayleigh numbers just above a critical value R^, we observed the forma- 
tion of chimneys, dendrite- free regions within the mushy- layer (Fig. |4jt). The instability 
controlling the transition from the no-chimney state to the mushy-layer mode of convec- 
tion depends on the mushy-layer Rayleigh number, dimensionlcss cell width A = VI/ k 
in addition to f <f, S and Ooo which depend on the liquid concentration Co- To identify 
a convective state in which chimneys form, we capture shadowgraph images to illus- 
trate two-dimensional fluid density variations for visualization of compositional plumes 
emerging from the mushy layer. We require a state of chimney convection to satisfy 
three criteria. First, there is a channel of zero solid fraction within the mushy region 
(Fig. 2k). Second, there is at least one vigorous compositional plume in the far- field 
liquid region. Since intermittent chimney p lumes (or convection in a breathing mode, 
[Sob mon fc Hartlevl[T99i IPeppin et aLll2008h are commonly observed in the experiment, 



we consider a chimney plume if its life time is longer than the diffusion time scale for 
boundary layer convection d 2 /D ~ 10 mins. Here d ~ 1mm is the average size of an 
isolated plume detached from the compositional boundary layer in the liquid region and 
D is the solute diffusivity. Third, other turbulent buoyant plumes in the near field next 
to the chimney plume have become substantially diminished which indicates suppressed 
boundary-layer mode convection in concert with the formation of the chimney. Transi- 
tion from a state of no flow to a state of chimney convection is readily identified through 
shadowgraph pictures. In Fig. [SJ we see the flow patterns in the liquid region when two 
sets of parameters (A,R m ) are chosen. Figure Hk shows that in a stable state of no flow 
(these data points are indicated by upward triangles in Fig. Wp)i there is no chimney flow 
in the far field of the liquid region in any of the cells, while small-scale turbulent plumes 
driven by the boundary-layer mode convection remain active. When the value of R m 
increases slightly (indicated by downward triangles in Fig. Wp), chimney flow penetrates 
into the far field in the two largest cells (Fig. ^jp) , which are readily discernible from the 
small-scale turbulent plumes and dominate the flow field in the liquid region. We thus 
conclude that the mushy-layers growing in the two largest cells are both in an unstable 
configuration that leads to a state of chimney convection. 

We determine the permeability of the mushy layer IIo in the following manner. By 
adjusting a single value of IIo that is the same for three different liquid concentrations, 
we determine values of R m from the experimental data so that the boundary between 
the stable and unstable regimes is closest to the stability curve derived from numerical 
calculations. As illustrated in Fig. [3 the resulting comparison with the theoretical curves 
is consistent with the observed stability to within one increment in cell width A. Hence, 
this comparison between theory and experiment allows us to empirically estimate the 
permeability as Ho = 1.6 xlO -5 cm 2 , independent of liquid concentration in the range 
24wt% < C < 27wt%. 

As seen in Fig. [7K, mushy-layer growth at a low initial liquid concentration (Co = 
24wt%) is mostly stable. Chimneys are present only at high values of R m and large A. 
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Figure 6. Shadowgraph images captured for Co = 25wt% and (a) 7=1.5K/min; (b) 
7=1.0K/min. Stability analysis results based on such images are presented in the regime di- 
agrams in Fig. 7b by downward and upward triangles, respectively. Both pictures are taken at 
Tb = — 40.0°C. Note the formation of a strong chimney plume in each of the two largest cells 
in (b), indicated by two white arrows. The horizontal white bars show 2 cm. 

Fig [7Jd shows that as the liquid concentration increases (to Cq = 25wt%), the critical 
R^, decreases. Moreover, the sensitivity of to A increases, abruptly increasing as the 
confinement increases (i.e. as A decreases) and preventing the onset of mushy-layer mode 
convection. When the liquid concentration further increases to Co = 27 wt%, solidification 
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Figure 7. Regime diagrams depicting the stability of the mushy-layer mode convection. Black 
symbols: mushy-layer grows in a stable state. Red symbols: unstable state with formation of 
at least one chimney channel within the mushy-layer, (a) Results for Co = 24wt%. (b) Results 
for Co = 25wt%. The triangles indicate the experimental results illustrated by shadowgraph 
pictures in Fig. [6] (c) Results for Co = 27wt%. In addition, similar experiments for Co = 23 wt% 
(not shown) resulted in no chimney formation for this range of cell widths and Rayleigh numbers. 
Experimental values of the Rayleigh number are calculated using the reference permeability 
IIo = 1.6x 10 -5 cm 2 for all three liquid concentrations, determined as described in the main text. 
Blue curves: numerical results for the stability curve indicating the critical RJjj as a function of 
A. 



of the sample with large A becomes unstable for all the values of R m obtained from the 
experiments as shown in Fig. [7fc. In summary, the system becomes less stable as the 
liquid concentration Co increases, which can be interpreted as follows. An increase in 
liquid concentration Co leads to the values of the concentration ratio ff, the Stefan 
number S, and the dimensionlcss liquid temperature 9^ all being reduced. Changes in 
these parameters influence the stab il ity of the system in different ways, as identified in the 
linear stability analysis of IWorsterl (Il992h . If is reduced, then the mean solid fraction 



is higher and hence the local permeability is reduced within the interior of the mush. In 
addition, a reduced value of leads to a thinner mushy layer and reduces the potential 
energy that is available to drive flow. Both of these effects make the system more stable. 
A smaller value of S implies that the influence of latent heat is smaller, so that the solid 
matrix is dissolved more easily. This promotes growth of the instability. A smaller value 
of #oo leads to a smaller heat flux from the liquid, so that the mushy layer grows thicker 
and provides greater potential energy to drive flow. The present theoretical simulations 
of finite amplitude convection reveal that the net combination of these competing effects 
results in the system becoming less stable as the liquid concentration Co is increased. 

4.3. Variability, multiple chimneys and morphological transitions between states 

Depending on the growth speed of the mushy-layer and the size of the sample-cell, several 
morphological transitions are observed in the states displaying chimney convection. To 
illustrate these morphological transitions we present these results as regime diagrams for 
three repetitions of experiments with Co = 27wt% in Fig. [5J The comparison between 
the three repetitions in panels (a), (b) and (c) is indicative of the role of fluctuations 
in influencing different experimental realisations. The stability boundary between states 
with and without flow is robust to within inter-experimental variation of one increment 
amongst the set of cell widths. The variability between observed states with and without 
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Figure 8. Results for the morphological states of the mushy-layer mode convection. Data ob- 
tained from three sets of experiments with Co = 27wt%. Black solid squares: stable state, 
mushy-layer with no-chimney; Red solid triangles: one chimney forms at the walls; Green 
open triangles: one chimney forms off the walls; Blue open circles: formation of two chimneys 
within the mushy-layer; Purple open squares: formation of three or more chimneys within the 
mushy-layer. Blue curves: numerical results for the stability curve indicating the critical R^ as 
a function of A (corresponding to a single convection cell with chimney at the wall.) 



flow under the same experimental conditions might be explained in part by the theo- 
retically suggested sub-critical instability behaviour: the form of perturbation required 
to trigger a transition to a finite-amplitude convecting state might only be achieved in 
a subset of the experiments. The variability may also be influenced by factors not de- 
scribed in the theory, such as weak heat loss through the walls or boundary-layer-mode 
convection. 

The regime diagrams in Fig. [8] also reveal a preference for chimney formation at the 
cell walls when the system undergoes a transition from a state of no flow to a state of 
chimney convection at A = A c i. Combining all the data in figure [8j for each cooling rate 
we determine the states at the smallest sample width that permits convection. Using this 
method we find four cell widths where a single chimney is observed at the cell wall, and 
the remaining two cell widths show either a chimney at the wall, or a chimney in the 
sample interior, for different repetitions under the same conditions. Localization of the 
mushy-layer mode of convection near the walls has be en observed in both experiments 
and models ( Giamei fc Kear 197fj Felicelli et al. 1998). Recent theoretical calculations 
show that chimney formation can be influenced by lateral heat transfer due to heat losses 
at the sidewalls. Sidewall heat transfer can have both direct th ermal effects, and indirec t 
effects by forming weak convective flows in the overlying fluid (jRoper et al.ll2007l . l201lh . 
The cooli ng effect of sidewal l heat losses can promote the formation of chimneys near to 
the walls (jRoper et a/.ll2007j) . Whilst localization of convection resulting from flow in the 
overlying fluid encourages chimney formation at the walls in three-dimensional mushy- 
layer cells, in two-dimensional flow the localisation promotes chimney formation away 
from the cell walls. Hence, this form of flow localisation cannot be the controlling factor 
in our quasi-two-dimensional Hclc-Shaw cell. As an alternative factor, our theoretical 
model suggests that the preference for wall chimneys close to the critical cell width for 
onset of convection may simply be a geometrical effect due to the finite-sample-width 
effects in our experiment. If the stability boundary A = A C) i for states with one chimney 
at the wall is interpreted as the minimum width of a single convection cell, then the 
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narrowest sample widths that can support convection necessarily have a chimney at the 
wall. A state with a central chimney separating two symmetric convection cells requires 
a minimum cell width width A = 2A c .i. Whilst we have not conducted detailed stability 
calculations for asymmetric states where an interior chimney is located away from the 
centre of the sample, it is plausible that certain forms of these states might be stable 
over the range A c i < A < 2A C i with the chimney separating a dominant convection cell 
of width greater than A Cj i from a narrower passive region. The possibility of asymmetric 
arrangements of convection cells is also likely to be of relevance in describing the observed 
transitions to states with two or more chimneys. For the same value of R m , transition 
from one-chimney states to two-chimney states occurred at a large spacing A Cj 2- Although 
the regime boundary for this transition varies slightly between different runs, we can 
conclude that A c ,2 > 2A c .i. A more detailed theoretical exploration of these morphological 
transitions presents an interesting opportunity for future work. 



5. Discussion 

In many practical contexts it is desirable to predict growth conditions which avoid 
the formation of chimneys in growing mushy layers. In cells that are sufficiently wide 
to permit the formation of multiple chimneys, experimental and theoretical models have 
been applied to infer a critical value R£, beyond which the mu shy-layer mode convection 
will occur (e.g. lFowlerlll985t I Worsted tl99a IChen et al.lll994l ). However, the cell width 
A pr ovides an ad d itiona l constraint on the existence of convective instability. In partic- 
ular, Wells et al. ( 2010l ) studied the evolution of finite amplitude states of convection 
with chimneys and showed that as A is decreased the chimneys are extinguished via a 
saddle- node bifurcation at some critical A c (R m ). The present work provides the first ex- 
perimental observation of this finite-samplc-size effect on buoyancy-driven convection in 
the solidification of mushy layers. 

We observe a discernible transition of the flow field in the liquid region overlying the 
mushy layer as mush-Raylcigh number R m and cell width are varied. The flow instability 
is systematically quantified by the sample spacing A and R m for a variety of liquid 
concentrations. In Fig. [6] we showed two sets of data from experiments conducted at 
two different cooling rates. At a high cooling rate, the resulting value of R m and the 
chimney spacing A (indicated by the upward triangles in figure [7b) characterize a stable, 
no-chimney state (Fig. [6^)- A small reduction of the cooling rate increases the value 
of R m (shown by the downward triangles in figure [7}d), and causes a flow instability 
in the mushy layer that results in a new state dominated by the formation of vigorous 
chimney plumes in the liquid region (Fig. |6}d). Using the same nominal experimental 
control parameters, the critical Raylcigh number in different runs varies less than 30%, 
as is also evident from the data shown in Fig. [8l Thus, these experimental studies provide 
reliable results and a general framework for direct comparisons with the prior numerical 
and theoretical studies, revealing that the system becomes less stable to the onset of flow 
if either the cell width, Rayleigh number or liquid concentration are increased. 

In determining the unknown reference permeability IIo of the mushy layer, we com- 
pared experimental and numerical results for the instability of the mushy-layer mode 
of convection summarized in Fig. [71 Good comparison can be achieved using a single 
value of the permeability IIo = 1.6xl0 -5 cm 2 across the entire experimental range of 
Co- This suggests that the reference permeability is relatively insensitive to changes in 
the liquid concentration. It is interesting to compare our estimate of the reference perme- 
ability in our quasi-two dimensional experiment to previous estimates. For examp le, for 
solidification in a fully three-dimensional experimental cell iTait fc Jaupartl (|1992h esti- 
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mated a mean permeability value of II = 8 x 10~ 4 cm 2 based on the Carmen-Kozeny 
relation and solid fractions inferred from measured temperature and salinity, whilst 
iNeufeld fc Wettlauferl (1200861) estimated a reference permeability TIq = 1.75 x 10 3 cm 



by com paring to a linear stab ility threshold for the onset of shear enhanced convection. 



Finally, I Chen &: Chenl (|1991I ) estimated a significantly lower value of the mean perme- 
ability II = 2.4 x 10~ 5 cm 2 based on the Carmen-Kozeny relation and solid fractions 
measured via computed tomography. However, this corresponds to a permeability both 
after the onset of convection and measured two days after the end of the experiment, 
allowing the system to evolve away from the initial value of permeability before con- 
vective onset. Hence, the difference in magnitude of the infe rred reference permeability 
betwee n our experiments and those of Tait fc Jaupartl (|l992l ) and INeufeld fc Wettlaufer 



(|20086l ) suggest that the permeability is significantly lower for solidification in a quasi 
two-dimensional cell than for a three-dimensional cell. 

Through this dynamical constraint of the permeability, the theoretical model for finite- 
amplitude convecting states successfully captures the form of the observed stability 
boundary to within the inter-expcrimental vari ability of one increm ent in sample width. 



This suggests that the theoretical framework of lWells ct al. (|2010l ) may provide a useful 



tool in future studies to investigate how chimney formation is influenced by additional 
factors, such as coupling to flow in the overlying liquid, in a fully nonlinear setting. 
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